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^ ■ ABSTRACT 

i Aims. We study the efficiency and reliability of cluster mass estimators that are based on the projected phase-space distribution of galaxies in a 
^> ■ cluster region. 

i Methods. We analyse a data-set of 62 clusters extracted from a concordance ACDM cosmological hydrodynamical simulation. We consider 
l/~) both dark matter (DM) particles and simulated galaxies as tracers of the clusters gravitational potential. Two cluster mass estimators are 
' considered: the virial mass estimator, corrected for the surface-pressure term, and a mass estimator (that we call M a ) based entirely on the 
| velocity dispersion estimate of the cluster. In order to simulate observations, galaxies (or DM particles) are first selected in cylinders of given 
radius (from 0.5 to 1.5 hr [ Mpc) and ^ 200 hT 1 Mpc length. Cluster members are then identified by applying a suitable interloper removal 
algorithm. 

Results. The virial mass estimator overestimates the true mass by a 10% on average, for sample sizes of ^ 60 cluster members. For similar 

■ sample sizes, underestimates the true mass by ^ 15%, on average. For smaller sample sizes, the bias of the virial mass estimator substantially 
Oh' increases, while the estimator becomes essentially unbiased. The dispersion of both mass estimates increases by a factor ~ 2 as the number 
^ of cluster members decreases from ~ 400 to ~ 20. 

5^ ■ It is possible to reduce the bias in the virial mass estimates either by removing clusters with significant evidence for subclustering or by selecting 
^ early-type galaxies, which substantially reduces the interloper contamination. Early-type galaxies cannot however be used to improve the M a 
C$ , estimates since their intrinsic velocity distribution is slightly biased relatively to that of the DM particles. 

' Radially-dependent incompleteness can drastically affect the virial mass estimates, but leaves the M a estimates almost unaffected. Other 
observational effects, like centering and velocity errors, and different observational apertures, have little effect on the mass estimates. 
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1. Introduction 

The masses of galaxy clusters are a very useful observable in 
cosmology. The number density of clusters of galaxies above a 
given mass threshold, and the evolution of this abundance with 
redshift z, provides in principle strong constraints on cosmo- 
logical models/parameters (see, e.g., Rosati et al. |2002 Voit 
2005 for recent reviews, and references therein). The advan- 
tage of using galaxy clusters is that they are massive, luminous 
objects, that can be detected relatively easily out to z ^ 1 by 
several techniques in different wavebands. The disadvantage is 
that they are rather complex objects, hence their masses are not 
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easily estimated, and can often be plagued by systematic effects 
uneasy to correct for. 

Since Zwicky's (1933 1 first estimate of a cluster mass, 
based on the application of the virial theorem to the projected 
phase-space distribution of galaxies in the Coma cluster, cluster 
mass determinations have always been taken with some cau- 
tion (see the historical review of Biviano 2000). A cluster mass 
estimate based on the observed projected phase-space distribu- 
tion of galaxies can be wrong because of several effects. Quite 
important, in this respect, are the projection effects leading to 
the inclusion of interlopers in the sample of presumed clus- 
ter members (see, e.g., Lucey 1983 Borgani et al. 119971 Cen 
1997, C97 hereafter). Biases in the cluster mass estimates can 
also occur when the studied cluster is far from virialization, 
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e.g. during the accretion phase of a massive group (see, e.g. 
Girardi & Biviano 2002 1, or if galaxies are biased tracers of the 
gravitational potential, which could happen as a consequence 
of dynamical friction (e.g. Biviano et al. 1 19921 Goto 2005), or 
as a consequence of infalling motions (e.g., Moss & Dickens 
119771 Biviano et al.[T997J. 

Because of the above mentioned problems, other methods 
of cluster mass determination have been considered, that are 
not based on the phase-space distribution of galaxies. While 
lensing mass estimates are also known to be affected by pro- 
jection effects, these are generally believed to be less severe 
(e.g. Reblinsky & Bartelmann IT999I RB99 hereafter; Clowe 
et al. 120041 but see Metzler et al. I1999> . unless the cluster is 
elongated or has substantial substructure along the line-of-sight 
(Athreya et al. 120021 Bartelmann & Steinmetz [l996| Gavazzi 
2005 , Oguri et al. 2005 1. Projection effects are even less impor- 
tant in the case of X-ray emission-based cluster mass estimates, 
since the X-ray emissivity is proportional to the square of the 
gas density. 

Problems with lensing and X-ray mass estimates do ex- 
ist, however. Masses determined with the lensing technique are 
affected by the mass-sheet degeneracy, that cannot always be 
broken (e.g. Dye et al. 2001 Cypriano et al. l2004> . The effect 
of intervening matter along the line-of-sight is to increase the 
weak-lensing mass estimates of clusters, especially at high red- 
shifts (Metzler et al. 1200 II Lombardi et al. 2005 , Wambsganss 
et al. 120041 12005 1. Lensing and X-ray-based mass estimates 
have often been found to be discordant, and this is generally in- 
terpreted as evidence of non-equilibrium (e.g. Wu 2000 , Clowe 
et al. l200Tl Athreya et al. l2002l Smith e t al. 120021 Cypriano et 
al. 120041 Ota et al. 120041 Bradac et al. 12005V In clusters un- 
dergoing merging events the X-ray luminosity and temperature 
can be boosted, thus leading to an overestimate of the cluster 
mass (e.g. Schindler & Muller [T993l Ricker & Sarazin I200T1 
Barrena et al. 2002 Diaferio et al. 2005 1. On the other hand, vi- 
olations of the condition of hydrostatic equilibrium, inaccurate 
modelling of the gas density profile and observational biases in 
the measure of the intra-cluster gas temperature, may lead to a 
sizeable underestimate of the cluster mass (e.g. Bartelmann & 
Steinmetz lT996l Kay et al. l2004l Rasia et al. l2004l 

A combination of several, independent cluster mass esti- 
mates is likely to provide the most accurate results. Moreover, 
with the Sloan Digital Sky Survey coming to completion (e.g. 
Abazajian et al. 2005 1, a large number of nearby clusters with 
> 30 galaxies with redshifts is now available (Miller et al. 
2005 1. For many of these clusters, X-ray data are not avail- 
able, and mass estimates must be based on optical data. Hence, 
mass estimates based on the projected phase-space distribution 
of galaxies are still very useful. 

Time is therefore mature for a re-assessment of the reliabil- 
ity of mass estimates of clusters based on the dynamics of their 
member galaxies. Previous analyses have generally considered 
only specific aspects of this topic. Initially, the reliability of 
different mass estimators has been assessed from N-body sim- 
ulations outside of a cosmological context (Danese et al. 119811 
Perea et al. 11990 Aceves & Perea 1999 1. In order to properly 
deal with this topic, clusters must however be identified within 
cosmological simulations with large box-sizes. This has been 



achieved by several studies (e.g. Frenk et al.Q990 Borgani et 
al. [T9971 C97; van H aarlem et al. [T997l hereafter vH97; RB99; 
Sanchis et al. 2004; Lokas et al. 20051, where, however, dark 
matter (DM hereafter) particles and not galaxies, were consid- 
ered as tracers of the potential. Van Kampen & Katgert ( 1997 > 
used the method of constrained random fields to increase the 
numerical resolution of their simulation, but by doing this they 
were unable to consider projection effects. In most analyses, in- 
terlopers were rejected using Yahil & Vidal's ( 1977 ) traditional 
3-cr-clipping method (but see, e.g., vH97 who also tested the 
more sophisticated method of den Hartog & Katgert 1996 and 
Sanchis et al. 2004 and Lokas et al. 2005 who tested their own 
methods of interloper removal). 

A general conclusion from these studies is that cluster mass 
estimates can be severely affected by projection effects. This 
happens mainly as a consequence of the cluster identifica- 
tion process in 2-d projected space, performed with Abell's 
( 119581 ) original algorithm. Frenk et al. ( 1990) argued that clus- 
ter masses are systematically over-estimated, but following 
studies (C97; vH97; RB99) concluded that cluster masses can 
be either over- or under-estimated, depending on the projection 
angle, the cluster mass, and the algorithm used to remove in- 
terlopers. Sanchis et al. (2004) and Lokas et al. (2005) found 
a rather good agreement between estimated and true cluster 
masses, when a rather large number of tracers of the potential 
was considered in each cluster. 

None of the above mentioned studies tried to identify galax- 
ies in the cosmological simulations. In those studies where this 
was achieved, it was generally found that the spatial distribu- 
tion of subhaloes, selected by their mass, is less concentrated 
than that of DM (e.g. Ghigna et al. [T998l Klypin et al. [T999l 
Gao et al. 2004). However, when gasdynamics was included in 
the simulations, galaxies, selected by luminosity, turned out to 
have a considerably more concentrated spatial distribution than 
subhaloes, and more similar to that of DM (Berlind et al. 2003 ; 
Gao et al. 2004, Nagai & Kravtsov 2005). This occurs because 
tidal stripping induce substantial mass loss from galaxy haloes, 
but very little stellar mass loss. 

As far as the velocity distribution of subhaloes is con- 
cerned, most studies have found it to be wider than that of DM 
particles (Colin et al. 2000; Diemand et al. 2004), at least near 
the cluster centre (Ghigna et al.|2000; Reed et al. l2005> . Based 
on hydrodynamical simulations, Frenk et al. (1996) however 
concluded that galaxies suffer significant dynamical friction, 
and are slowed down relatively to DM particles, to such an ex- 
tent that cluster mass estimates based on the velocity dispersion 
of galaxies are likely to be in error by factors of 0.25-0.75, de- 
pending on the masses of the galaxies selected as tracers. On 
the other hand, Berlind et al. (2003) found that galaxies have 
only a mild velocity bias with respect to DM, and Faltenbacher 
et al. (2005) found that galaxies, if anything, move slightly 
faster than DM particles. 

In this paper the reliability of cluster mass estimates based 
on the dynamics of their member galaxies is reconsidered on 
the basis of a set of clusters extracted from a large cosmo- 
logical hydrodynamical simulation (Borgani et al. 2004| see 
also Murante et al. 12004b . performed using the TREE+SPH 
GADGET-2 code (Springel et al. l200ll Springel 120051 . for a 
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concordance ACDM model. This simulation samples a fairly 
large volume (box size of 192 1 Mpc, where h is the Hubble 
constant in units of 100 km s _1 Mpc -1 ) which contains more 
than 100 clusters with mass above 10 u h~ l M Q , while also hav- 
ing good enough resolution (todm = 4.62 x 10 9 h M® for the 
mass of the DM particles) to allow resolving halos hosting 
bright galaxies. Finally, the inclusion of the processes of ra- 
diative cooling, star formation and supernova feedback allows 
us to have a realistic description of the gas evolution and of the 
galaxy formation process. 

In this paper we address the question of how accurate are 
cluster mass determinations based on the dynamics of their 
member galaxies, under a variety of observational conditions, 
and for a wide range of cluster masses, but independently of the 
clusters identification procedure. No attempt is made here to 
simulate the observational identification algorithms of galaxy 
clusters, such as Abell's ( 1958 1 original one. In this sense, our 
approach is different from previous ones (Frenk et al. 119901 
C97; vH97; RB99), in that we disentangle the problem of clus- 
ters mass estimation from that of clusters identification. The ra- 
tionale for this choice is that today there is no standard cluster 
identification algorithm. Automated scans of digitized plates 
have since replaced Abell's eyeball identification of galaxy 
clusters (Dalton et al. 1 19921 Lumsden et al. 1 19921 Lopes et al. 
2004 1, and much more sophisticated algorithms than Abell's 
( 1958 1 have been applied for extracting the 2-dimensional sig- 
nal produced by a galaxy overdensity (e.g. Ramella et al. 120011 
Gladders & Yee 2005 and references therein). Moreover, those 
clusters whose masses are derived using member galaxies, are 
not necessarily optically selected (e.g. Popesso et al. 2005 1. 
The results of this paper are therefore useful for a better un- 
derstanding of why different techniques (lensing, X-ray, galax- 
ies) can lead to discrepant cluster mass estimates, and also for 
the study of scaling relations of different cluster properties with 
cluster masses. Translating the results of this paper to the study 
of the distribution of observed cluster masses in a given survey, 
requires convolution with the selection function of the survey 
itself. 

The plan of this paper is as follows. In Sect.|2]the set of sim- 
ulated clusters and its characteristics are described. In Sect. [5] 
the steps involved in the determination of cluster dynamical 
quantities are briefly described. The analyses of the cluster 
dynamics in projected phase-space are described in Sect. 0] 
Results are discussed and summarized in Sect. [5] In Sect. [6] 
we provide our conclusions. 

2. The simulated clusters 

The set of simulated clusters, analysed in this paper, are ex- 
tracted from a large cosmological hydrodynamical simulation 
presented by Borgani et al. (2004i, and performed using the 
massively parallel Tree+SPH GADGET-2 code (Springel et al. 
I200TI Springel 120051 . We refer to Borgani et al. (2004i for 
a complete description of this simulation and of the cluster 
identification procedure. We provide here their main charac- 
teristics and describe in some details the galaxy identifica- 
tion procedure. The simulation assumes a cosmological model 
with O = 0.3, f2 A = 0.7, Q b = 0.019 h~ 2 , h = 0.7, and 



o- 8 = 0.8. The box size is 192ft" 1 Mpc. We used 480 3 DM 
particles and (initially) as many gas particles, for a mass reso- 
lution of moM = 4.62 x 10 9 /j _1 M g . The Plummer-equivalent 
softening length was set to e = 7.5 h~ l kpc at z — 0. The simu- 
lation code includes explicit energy end entropy conservation, 
radiative cooling, a uniform time-dependent UV background 
(Haardt & Madau 1996 1, the self -regulated hybrid multi-phase 
model for star formation (Springel & Hernquist 2003 1, and a 
phenomenological model for galactic winds powered by Type- 
II supernovae. 

We identify galaxy clusters at redshift z — by applying a 
standard Friends-of-friends (FoF) analysis to the DM particle 
set, with linking length A = 0.15 in units of the mean interpar- 
ticle distance. After the FoF identification, we applied a spheri- 
cal overdensity criterion to find the virial region of each cluster 
(corresponding to an overdensity of « 100 times the critical 
density for the adopted cosmology). We identify 117 galaxy 
clusters whose virial mass is larger than 10 14 h~ l M Q within our 
simulated volume. A subset of these clusters is selected for the 
analysis discussed in the following. The position of the mini- 
mum potential particle of each FoF group is then used as the 
centre for a spherical overdensity algorithm, which identifies 
the radii encompassing different overdensities. 

Galaxies are identified using the publicly available algo- 
rithm 1 SKID (StadelEUnD- SKID calculates a density field us- 
ing all the particles. Density is assigned with a spline-kernel 
similar to that employed by SPH codes; an important parameter 
is the number of neighbouring particles N sm on which the den- 
sity is softened. Particles are then associated to the local max- 
ima of such field, which should represent the positions of the 
substructure. To this aim, particles are "moved" along the gra- 
dient of the field until they begin to oscillate around the peaks 
and are then grouped using a FoF analysis on the "moved" po- 
sitions. An unbinding procedure is finally applied on the result- 
ing FoF groups, to discard particles which are not gravitation- 
ally bound to the group. All particles within a sphere of given 
radius are used to evaluate the gravitational potential. This step 
is performed on the "true" (i.e. not "moved") positions. All the 
length scales involved in the process are set starting from the 
parameter t which represents the typical size of the objects to 
be found. Apart from the evaluation of the density field, it is 
possible to select to which types of particles (DM, stars, and 
gas) the algorithm must be applied. We used only star particles 
to identify our galaxies. 

The problem of reliably detecting sub-structures within 
given DM haloes is long-lasting and still not uniquely solved. 
It is outside the purpose of the present work to discuss such an 
issue. We performed a number of tests on the SKID galaxies, 
to verify that the least possible number of objects is missed by 
the analysis and that the identified object are real ones. Such 
tests are presented elsewhere (Murante et al, in preparation). 
Our galaxies have been identified by using t = 20/z~' kpc, 
which, after a trial-and-error procedure, turned out to be the 
optimal choice. This choice is also motivated by the fact that 
20/i~' kpc roughly correspond to the effective force resolution 

1 See 

http : //www-hpcc . astro . Washington . edu/tool s/skid . html 
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Fig. 1. The projected number density profiles of simulated 
galaxies (solid lines) and real galaxies from the ENACS sam- 
ple (dots with error bars). The three solid lines in each plot 
correspond to three orthogonal projections. The normalisations 
of the simulated galaxies profiles have been arbitrarily scaled 
to match the observed profiles. The clustercentric distances on 
the x-axis are in units of the cluster virial radius. Top panel: 
simulated galaxies with formation redshift zj > 1.25 (see text) 
vs. early-type galaxies from the ENACS sample. Bottom panel: 
simulated galaxies with zj < 1.25 vs. late-type galaxies from 
the ENACS sample. 
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Fig. 2. The line-of-sight velocity dispersion profiles (in units 
of the total cluster velocity dispersion) of simulated galaxies 
(solid lines) and real galaxies from the ENACS sample (dots 
with error bars). The clustercentric distances on the x-axis are 
in units of the cluster virial radius. The three solid lines in each 
plot correspond to three orthogonal projections. Top panel: 
simulated galaxies with formation redshift z/ > 1.25 (see text) 
vs. early-type galaxies from the ENACS sample. Bottom panel: 
simulated galaxies with zj < 1.25 vs. late-type galaxies from 
the ENACS sample. 



of the simulation. We also determined that distinct SKID anal- 
yses are needed, with different values of N sm , that we assumed 
to be 16, 32 and 64. Some galaxies are "missed" using only 
one value for N sm . Therefore, we built our catalogue by joining 
together the results of the three SKID analyses, with the fol- 
lowing rule: if a star particle is a member of a galaxy for one 
value of N sm , it is considered to belong to that galaxy; if a parti- 
cle is member of two different galaxies for two different values 
of N sm , the two galaxies are then merged into a single object. 
SKID objects with less than 32 particles are discarded. 

For the purpose of this paper it is not useful to consider 
clusters with less than 10 galaxies, hence in the following we 
restrict our analysis to those 62 clusters with at least 10 galaxies 
within a sphere of 'virial' radius r v , defined as the radius where 
the mass density of the cluster equals 200 times the critical 
density of the Universe. 

In Borgani et al. (2004i we have shown that these simulated 
clusters have on average a star fraction which is ~ 50 per cent 
larger than the typical observed values (Lin, Mohr & Stanford 
2003 1. However, this overcooling is mainly contributed by the 
central galaxy. Therefore, our simulations are expected to reli- 
ably describe the way in which cluster galaxies trace the under- 
lying cluster dynamics. 

For each galaxy, we define its formation redshift, tf, as the 
average value of the formation redshifts of all its member star 
particles. The value of Zf is also used to classify our galax- 
ies into "early-type" and "late-type". Specifically, we consider 
a galaxy to be of early-type if z/ > 1.25, and of late-type if 
Zf < 1.25. This choice is admittedly rather crude and is based 
on the comparison of our simulated galaxy sample with the 



sample of cluster galaxies of the ESO Nearby Abell Clusters 
Survey (ENACS, Katgert et al. IT9961 [1998). In the ENACS, 
64% of the galaxies identified as cluster members within 1 .5 r v 
are classified as early-type galaxies (ellipticals, lenticulars, or 
an intermediate class between these two; see Biviano et al. 
2002 Thomas & Katgert 2006 1. Similarly, in the simulated 
clusters, 64% of the galaxies within 1.5 r v have zj > 1.25. 

Further support to the identification of Zf > 1-25 galaxies 
as early-type galaxies comes from the comparison of their dis- 
tribution in projected phase-space with that observed for the 
ENACS cluster galaxies. Since the number of galaxies per sim- 
ulated cluster is rather limited, clusters are stacked together, by 
scaling the galaxy clustercentric distances and velocities by the 
value of r v and the velocity dispersion of the cluster they be- 
long to. A similar procedure is applied to the observed galaxies 
(see Katgert et al. 120041 and Biviano & Katgert 2004). The pro- 
jected number density profiles, ~E{R), and the line-of-sight ve- 
locity dispersion profiles, <t p (R), of real and simulated galaxies 
are shown in Figures^andlJl respectively, separately for early- 
and late-type galaxies (top and bottom panels, respectively). 
The simulated galaxies £(/?)'s have been arbitrarily rescaled to 
match the observed £(/?)'s, since here we are only interested in 
the relative distributions of simulated and real galaxies. 

The profiles of observed and simulated galaxies are rather 
similar, both for the early- and the late-types. The main differ- 
ence between the observed and simulated profiles is the lower 
number density and smaller velocity dispersion of the simu- 
lated early-type galaxies in the inner regions (R/r v ^ 0.1). This 
difference is caused to the well-known overmerging problem 
which affects the simulations in the denser regions. While this 
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can be a source of concern when comparing simulated vs. ob- 
served cluster mass profiles, we argue that the analysis of clus- 
ter masses is not significantly influenced by this problem. In 
fact, the apparent underdensity of simulated early-type galaxies 
in the inner regions, as compared to the real ones, is not really 
significant. The fraction of early-type galaxies within the virial 
region, that lie within 0.1 r v is 12.9% for the simulated galax- 
ies, and 14.5% for the observed ones. As a consequence, no 
significant effect is expected on cluster harmonic mean radius 
estimates (which enters the virial mass estimate, see Sect. 0, 
also because these fractions are low. For the same reason, no 
significant effect is expected on cluster global velocity disper- 
sion estimates, despite the fact that the difference in the veloc- 
ity dispersions of simulated and observed early-type galaxies is 
significant within 0.1 r v . In fact, the ratio between the velocity 
dispersion estimates obtained using all galaxies and only those 
with R > 0.1 r v is 0.994 for the simulated galaxies, and 1.008 
for the real ones. 

We conclude that the similarity between the profiles of sim- 
ulated and observed galaxies lends support to our choice of 
Zf = 1.25 for separating early- and late-type galaxies in our 
simulations. 




1000 

(7 v (km/s) 



Fig. 3. The true mass vs. the 3-d velocity dispersion (both com- 
puted within a sphere of radius r v ) for the 62 simulated clusters 
with at least 10 galaxies within r v . The best fitting cubic law 
relation is shown (solid line). 
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3. The cluster mass estimates 

Several different definitions of 'virial mass' have been given in 
the literature. This term may be used to define the total mass 
of a cluster within a radius of given overdensity, i.e. the 'virial 
radius', typically the radius where the mean cluster mass den- 
sity equals 100 or 200 times the Universe critical density (see, 
e.g., Lokas & Mamon 2003 1. On the other hand, the same term 
is used to define the mass estimated through application of the 
virial theorem to the cluster galaxies within an observationally 
defined aperture (e.g. Biviano et al.Q993; Girardi et al. 1998 1. 
This estimate requires not only an estimate of the cluster ve- 
locity dispersion, but also of the harmonic mean radius of the 
spatial distribution of cluster galaxies. To further complicate 
the issue, twice the harmonic mean radius is also usually re- 
ferred to as the 'virial radius' (e.g. Girardi et al. 1998 , Merchan 



&Zandivarez 2005). 

The virial mass estimate requires correction for the surface 



pressure term (The & White 1986), unless the entire system 
is contained within the observationally defined aperture radius 
(see also Maccio et al. 2003 ). Unfortunately, it is still relatively 
uncommon to see this correction applied (see, e.g., Koranyi et 
al.[T998l. Neglecting this correction leads to overestimate the 
mass of a system (see, e.g., Carlberg et al. 1997ai, and this can 
partly account for some of the claimed discrepancies between 
optically and X-ray derived cluster mass estimates. 

For the sake of clarity, we detail in the following the proce- 
dure that we apply to our simulated clusters in projection in or- 
der to simulate an observational estimate of their masses. Such 
a procedure has recently been applied to a large set of nearby 
clusters by Popesso et al. ( 20051 12006V 

We define a cluster 'true' mass, M v , as the total mass within 
the radius r v , where the mass density of the cluster equals 200 
times the critical density of the Universe at the cluster redshift 
(we call cr v the 3-d velocity dispersion of the DM particles 
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Fig. 4. The average integrated velocity dispersion profile for 
the 62 simulated clusters. The 3-d velocity dispersion within a 
sphere of radius a is shown as a function of a/r v , in units of cr v . 
1-cr error bars are shown. 

within the same radius). The virial mass estimate M v within 
the same radius is computed as follows. 

1 . Define an aperture radius, a, within which to perform the 
dynamical analysis, and select all DM particles or all galax- 
ies within a cylinder of radius a and ^ 200 h Mpc length, 
along each of three orthogonal projections. Cluster centres 
are computed as the positions of the minimum potential 
particle in each cluster (see Sect. |2}, unless otherwise indi- 
cated (see Sect. 14. 5> . In virtually all cases, these centres cor- 
respond (to within ~ 10 kpc) with the peaks of the cluster 
X-ray emissivity. The same centre choice is usually adopted 
observationally, when X-ray data are available (see, e.g., 
Biviano et al. 119971 Popesso et al. l2005> . 

2. Select cluster members. This is done first by using a cut in 
line-of-sight velocity space of ±4000 km s _1 with respect to 
the mean cluster velocity, initially defined by applying the 
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biweight estimator 2 (see Beers et al. 1990 1 to all the objects 
within 0.5 hr x Mpc from the cluster centre. In a second 
step, the weighted gap procedure of Girardi et al. ( 1993 1 is 
applied on the remaining objects. On the objects that pass 
the weighted gap selection we finally apply Katgert et al.'s 
( 2004 1 procedure that makes use of the location of galaxies 
(or particles) in projected phase-space (see also den Hartog 
& Katgert's[l99lj). 

Determine the 'robust' estimate of the line-of-sight velocity 
dispersion cr a p , within the aperture a, using the biweight or 
gapper estimator, depending on the number of data avail- 
able, > 15 or < 15, respectively (see Beers et al. 1990 1. 
Correct the velocity dispersion estimate for velocity errors 
(if these are added to the simulated data, see Sect. 14. 5> fol- 
lowing the prescriptions of Danese et al. ([1980 1. 
Determine the projected harmonic mean radius r a4 ,, within 
the aperture a. 

Obtain a first estimate of the mass from 



M a = 3ncrl^ p /G, 



(1) 



where G is the gravitational constant, the factor 3^/2 is the 
deprojection factor, see Limber & Mathews (1960 1, and a 
factor 2 is needed to convert the harmonic mean radius into 
the 'virial radius' of Girardi et al. ( 1998 1. 



Estimate the Navarro et al. ( 1997 NFW hereafter) concen- 
tration parameter using the relation c = 4(<x ai/J /700)~ ' 306 . 
The normalization of the relation is taken from Katgert et 
al. i 



(|2p04 1, and the exponent of the relation is derived from 
Dolag et al. (|2004l, under the assumption that the mass 
scales with the third power of the velocity dispersion. 
Correct the mass estimate for the surface pressure term, f sp , 
obtained for a NFW profile with concentration c: M fl c = 
f sp M„, assuming isotropic orbits (see eq.(8) in Girardi et 
al. 1998 1. Note that the value of the correction factor is not 
negligible (f sp — 0.84 on average for our simulated clusters 
observed out to an aperture radius a = 1.5 h Mpc). 
Determine r v , an estimate of r v , as f v = a[p a /(2Q0p c (z))] 1 ^, 
where p c (z) is the critical density of the Universe at the clus- 
ter redshift, p a = M a , c /(4na 3 /3), and £ is the local slope of 
a NFW profile of concentration c at the radius a. 
Determine M v , an estimate of M v , by extrapolating or inter- 
polating M a>c from a to r v using a NFW profile with con- 
centration c. 



Quite often, cosmological constraints have been obtained 
directly from the distribution of cluster velocity dispersions, 
used as proxies for the cluster masses (see, e.g., Girardi et al. 
I1993> . It is therefore useful to also consider an alternative (and 
simpler) mass estimator entirely based on the line-of-sight es- 
timate of the velocity dispersion within a given aperture, cr a p . 
The relation between M v and cr v for our simulated clusters is 
shown in Fig. [3] The line in that plot represents the best-fitting 
cubic relation 



MvS W_^_f 

UO^kms- 1 / 



x 10 14 /r'Mp, 



(2) 
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2 Here and throughout this paper we use the biweight estimator for 
the average and dispersion, unless otherwise indicated. 
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Fig. 5. Verifying the performance of the mass estimators: the 
comparison between the masses estimated using full phase- 
space information, and the true cluster masses. Top-left panel: 
ratio between the virial and true mass, M v /M v vs. the true mass, 
M v . Bottom-left panel: ratio between the mass obtained from 
cr v and the true mass, M a /M v vs. the true mass M v , Top-right 
panel: distribution of M v /M v , Bottom-right panel: distribution 
of Ma-IMy. The masses M v and Mg- are estimated from the 
full phase-space distributions of 5000 DM particles per sys- 
tem (when available; otherwise all DM particles are selected) 
within the cluster virial radius. 

where A = 1.50 + 0.02 (1-cr error). We note that this rela- 
tion is approximately independent of the cosmological model 
(Borgani et al. 1999 1. It is possible to use this relation to obtain 
a mass estimate entirely based on the velocity dispersion esti- 
mate. Since can differ significantly from <x,, depending 
on the aperture radius a chosen by the observer, we need to 
apply a correction that depends on the ratio a/r v . For each of 
the 62 clusters, we determine the 3-d velocity dispersion cr a , 
within a sphere of radius a, for different values of a. The aver- 
age cr a /cr v vs. a/r v profile of our 62 clusters is shown in Fig.|4] 

Hence, given a line-of-sight velocity dispersion estimate 
within a given observational aperture, <x„ p , we proceed as fol- 
lows. We use y3cr a ,p to determine an initial estimate of the 
mass of the system through eq. |2] We then obtain an estimate 
f v of r v , following steps 7 and 9 above (note that step 8 is not 
needed in this case, because the relation provided by eq. pre- 
lates the velocity dispersion to the true cluster mass, for which 
no surface term correction is needed). Using the cr a /cr v vs. a/r v 
profile of Fig. 0] and replacing the true quantities r v and cr a , 
with their estimates r v and y3cr a ,p, respectively, we obtain an 
estimate of cr v which we use to determine the M a estimate of 
M v through eq.|2 

In order to verify the performance of these two mass esti- 
mators, we have applied them to the clusters in the simulation, 
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making use of the full phase-space information. In this case, 
there is no need to apply the interloper-rejection procedure out- 
lined above (step 2), since we randomly select 5000 DM par- 
ticles in each cluster within a sphere (not a cylinder) of radius 
a = r v (when this sphere contains less than 5000 DM particles, 
all particles are selected). Also the de-projection factor 37r/2 
(step 6) is not needed. As for the estimate of M a , we directly 
apply eq.|2l In Fig.[5]we show the ratio between the virial (M v ) 
and the true masses (upper panel), and M, T and the true masses 
(lower panel) for our sample of 62 clusters. On average we find 
M v /M v = 1.03 with a dispersion of 0.09, and M^/M, = 1.00 
with a dispersion of 0. 12 (see also Tormen et al. 1997 for a sim- 
ilar analysis). Perfect identity is not expected, since clusters are 
not fully virialized, but the good agreement indicates that the 
estimators are unbiased and the clusters in the simulation are 
close to virialization within r v , as expected for z ~ clusters in 
a ACDM cosmology. 

4. Dynamical analysis in projected phase-space 

In order to mimic observations, the simulated clusters are 
placed at a distance of 250 h~ x Mpc from the hypothetical ob- 
server. Three orthogonal projections are considered for each 
cluster, leading to a total of 186 (62 x 3) cluster projections 
(except when we remove the cases with significant evidence 
for subclustering, see Sect. 14. 2\ . DM particles or galaxies are 
then selected within cylinders of given radius, chosen to be a 
fraction (from 1/3 to 1) of an Abell radius. This is meant to 
reproduce typical observational procedures where the virial ra- 
dius of the system is unknown prior to observation. The length 
of the cylinder is set by the simulation cube. 

In the following we analyse how cluster mass estimates 
are affected by several observational effects. In particular, we 
consider the effects of different sample sizes, different obser- 
vational apertures, observational errors, incompleteness, and 
subclusters detected in projected phase-space. In principle, it 
would be desirable to perform all these analyses using the sim- 
ulated galaxies as tracers of the potential, rather than the DM 
particles. As a matter of fact, the analyses performed on the 
simulated galaxies are likely to provide a reliable picture of the 
real observational situation, since the phase-space distributions 
of simulated galaxies are not too different from the observed 
ones (see Sect. |3J. On the other hand, the DM particles are 
not distributed like the galaxies in the simulated clusters (see 
Sect. 14. 61 . Unfortunately, however, there are only few galaxies 
per cluster in our simulated set, hence small-number statistical 
noise becomes a dominant source of scatter. In order to over- 
come this problem, we have also considered DM particles as 
tracers of the gravitational potential, even if, of course, these 
are not true 'observables'. 

The most relevant results of the dynamical analyses in pro- 
jected phase-space are summarized in Table ^ The different 
analyses performed are identified by their number in the 1st 
Col. of Tabled Col. 2 lists the observational aperture size, a, 
in hT x Mpc, i.e. the radius of the cylinder within which particles 
or galaxies are selected. Col. 3 lists the number of DM particles 
per projection initially selected for the analysis within a cylin- 
der of given radius. In the case of galaxies, all galaxies (of the 



specified type) are always selected, and this number changes 
from projection to projection, hence we list the average over 
all projections. This is also the case for DM particles, when 
we select an identical number of DM particles as of available 
galaxies within the cylinder (see Sect. 14. 6> . In Col. 4 'DM' or 
'G' indicates whether DM particles or Galaxies are chosen as 
tracers of the potential. When only galaxies of a given type are 
selected, this is further specified after the letter 'G'. Col. 5 lists 
the type of selection performed on DM particles, which can 
either be random ('Rdm') or of a type that reproduces a single- 
mask, multi-slit observation ('Slit'). In the case of galaxies, 
'All' indicates that all galaxies in the cylinder are selected. The 
note 'projections with subclusters excluded' indicates that the 
definition of the sample used is the same as in the line above, 
except that cluster projections with significant evidence of sub- 
structure have been removed from the sample. Col. 6 and 7 
list the average of the number of DM particles or galaxies se- 
lected as cluster members within the chosen aperture, and the 
fraction f OM of these that actually lie outside the sphere of ra- 
dius equal to the chosen aperture, and that we call "interlopers". 
Cols. 8 and 9 list the average ratios of the estimated-to-true ve- 
locity dispersion V3<x fl _ p /cr a , and of the estimated-to-true har- 
monic mean radius |r a / ,/r fl , respectively, all quantities being 
estimated within the given aperture a (note that the projected 
quantities are multiplied by the deprojection factors). Col. 10 
lists the average ratio of the estimated-to-true virial radius, 
f v /r v . Cols. 11 and 12 list the average ratios of the two mass 
estimates (the virial estimate and the velocity-dispersion based 
estimate, respectively) to the true mass, M v /M v and Mg-/M v , 
respectively. The dispersions of the quantities listed in Cols. 
6-12 are given (within brackets) in the line immediately be- 
low. Cols. 13 and 14 list the fractions of cluster projections for 
which the virial mass is a strong under- or over-estimate of the 
true mass, MJM V < 1 /2, and M v /M v > 2, respectively. 

4.1. Projection effects 

In order to quantify the importance of projection effects on 
cluster mass estimates, we first consider a large number (500) 
of reliable tracers of the gravitational potential (DM particles), 
with no errors on their velocities, positions, and no radial in- 
completeness bias (see Tabled Id- no.l). The 500 DM particles 
are randomly selected along each of three orthogonal line-of- 
sights for each cluster, within a cylinder of 1.5 hr x Mpc radius. 
Interlopers are then rejected using the procedure described in 
Sect. |3 The case studied here is of course an idealized situa- 
tion never achieved observationally, but this analysis serves to 
emphasize the importance of projection effects on cluster mass 
estimates, independently of any other observational bias. 

On average, both the virial mass estimate and the M lT esti- 
mate are accurate to ^ 10%, but the scatter around the true mass 
values for the different projections is quite large (see Fig. [6). 
The virial estimate somewhat overestimates the true mass, and 
this is caused by an overestimate of the harmonic mean ra- 
dius, while Mu- somewhat underestimates the true mass, and 
this is caused by an underestimate of the velocity dispersion 
(see Table[0. 
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Table 1. The ratios between the estimated dynamical quantities in projection and the true ones. Quantities in brackets are the 
dispersions on the estimated quantities. 



Id. 


a 

(hr x Mpc) 


N 


Tracer 


Sel. 


N m 


foul 


^0~a,p/0- a 


2 ^"a,p 1 


r v /r v 


M r /M r 


MJM r 


/<0.5 


f>2 


1 


1.50 


500 


DM 


Rdm 


402 


0.18 


0.93 


1.22 


1.08 


1.10 


0.88 


0.01 


0.03 












(48) 


(0.11) 


(0.10) 


(0.12) 


(0.10) 


(0.31) 


(0.28) 






2 


projections with subclusters excluded 


396 


0.17 


0.90 


1.21 


1.06 


1.02 


0.81 


0.00 


0.00 












(50) 


(0.10) 


(0.10) 


(0.11) 


(0.08) 


(0.25) 


(0.24) 






3 


1.50 


100 


DM 


Rdm 


80 


0.18 


0.93 


1.20 


1.07 


1.10 


0.81 


0.04 


0.10 












(9) 


(0.11) 


(0.14) 


(0.19) 


(0.14) 


(0.46) 


(0.38) 






4 


1.50 


20 


DM 


Rdm 


17 


0.20 


0.89 


1.37 


1.11 


1.25 


0.83 


0.09 


0.25 












(2) 


(0.14) 


(0.27) 


(0.28) 


(0.25) 


(0.79) 


(0.58) 






5 


1.50 


500 


DM 


Slit 


198 


0.26 


0.90 


1.63 


1.18 


1.51 


0.83 


0.00 


0.16 












(51) 


(0.13) 


(0.12) 


(0.17) 


(0.10) 


(0.41) 


(0.28) 






6 


0.50 


500 


DM 


Rdm 


459 


0.28 


0.94 


1.20 


0.97 


1.08 


0.87 


0.01 


0.02 












(46) 


(0.09) 


(0.10) 


(0.07) 


(0.09) 


(0.32) 


(0.29) 






7 


1.50 


39 


G 


All 


30 


0.26 


0.95 


1.39 


1.16 


1.40 


0.94 


0.11 


0.27 












(12) 


(0.15) 


(0.27) 


(0.25) 


(0.24) 


(0.85) 


(0.62) 






8 


1.50 


39 


DM 


Rdm 


29 


0.24 


0.98 


1.37 


1.17 


1.46 


1.03 


0.03 


0.31 












(12) 


(0.14) 


(0.26) 


(0.26) 


(0.24) 


(0.89) 


(0.68) 






9 


1.50 


19 


Gz f > 1.25 


All 


16 


0.18 


0.89 


1.21 


1.04 


1.01 


0.80 


0.20 


0.15 












(8) 


(0.15) 


(0.30) 


(0.30) 


(0.26) 


(0.74) 


(0.64) 






10 


1.50 


19 


DM 


Rdm 


15 


0.25 


0.93 


1.41 


1.13 


1.32 


0.92 


0.14 


0.34 












(7) 


(0.18) 


(0.37) 


(0.32) 


(0.35) 


(1.14) 


(0.83) 







In order to better understand the effects of projection, we 
examine the projected phase-space distribution of a stacked 
cluster in some detail. The stacked cluster is built as follows. 
In each cluster of the simulation set we randomly select 500 
DM particles along three orthogonal projections, and identify 
the cluster members. We then normalize the clustercentric dis- 
tances of these particles by the virial radius of the cluster they 
belong to, and their velocities (relative to their cluster mean ve- 
locity) by the velocity dispersion of the cluster they belong to, 
computed within one virial radius. We then stack the 62 clusters 
using the normalized radii and velocities. Finally, we return to 
physical units by multiplying the normalized radii and veloci- 
ties of the DM particles in the stacked cluster, respectively by 
the average virial radius (0.99 h Mpc), and by the average 
velocity dispersion (657 km s~') of the 62 clusters. 

In Fig.0we show the projected phase-space distribution of 
a randomly selected subset of the DM particles in the stacked 
cluster. Three orthogonal projections are stacked together in the 
Figure. Different symbols are used to represent: i) DM particles 
contained within a sphere of 1.5 h~ x Mpc radius and selected 
as cluster members by the interloper rejection procedure (filled 
squares), ii) DM particles located within this same sphere yet 
incorrectly rejected as interlopers (open squares), and iii) DM 
particles located outside this same sphere and yet selected as 
cluster members (x symbols). Clearly, the interloper rejection 
technique excludes only a few DM particles that are in fact 
within the 1.5 hr l Mpc sphere, but keeps several DM particles 
that are outside the 1.5 h~ l Mpc sphere but have projected ve- 
locities that make them indistinguishable from the real cluster 
members in the projected phase-space diagram (what we have 
called interlopers, see Sect.|4}. 

The underestimate of the velocity dispersion of the cluster 
is partly due to the rejection of real cluster members at rel- 



atively high velocities. Most of the underestimate is however 
due to the inclusion of interlopers that are currently infalling 
toward the cluster along a filament. These interlopers lie within 
£ 10 h l Mpc from the cluster centre (more distant DM par- 
ticles have quite a different recession velocity from the cluster 
mean and are rejected by the interloper technique) and their 
infall motion makes their velocities resemble the average clus- 
ter velocity. In addition, they have an intrinsically small ve- 
locity dispersion within their filament. As a consequence, the 
projected velocity dispersion of these interlopers turn out to 
be smaller than that of the cluster within the selected aperture. 
On the other hand, their spatial distribution is not as centrally 
concentrated as that of the cluster members within the selected 
aperture, as shown in Fig. [S] This causes the harmonic mean 
radius to be overestimated. 

Therefore, interlopers cause the overestimate of the har- 
monic mean radius, and, at the same time, the underestimate of 
the velocity dispersion, as already noted by C97 and Diaferio et 
al. ([1999 ). Since the strongest of the two effects is the harmonic 
mean radius overestimate, also the virial mass is overestimated. 
On the other hand, provides an underestimate of the mass 
because it only depends on the velocity dispersion estimate. 

Despite these problems, the interloper selection technique 
seems to work reasonably well. On average only ~ 1% of the 
real cluster members are rejected as interlopers, while on av- 
erage 18% of the selected members are unrecognized interlop- 
ers (see f om in Table [Q. However, projection effects can be- 
come critical for individual cases. In 3% of the cases the clus- 
ter masses are overestimated by a factor > 2, and in 1% of 
the cases they are underestimated by a factor < 1/2, using the 
virial mass estimator (see Fig. [6] and the last two columns in 
Table Large errors in the virial mass estimates correspond 
to extreme cases of failure of the interloper rejection algorithm, 
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Fig. 6. The comparison between the masses estimated using 
only projected phase-space information and the true cluster 
masses. Top-left panel: the ratio between the virial and the true 
mass, M v /M v vs. the true mass M v . Filled dots identify those 
projections where no significant evidence for substructure is 
found at the 0.01 c.l. Top-right panel: distributions of M v /M v 
for the projections with and without significant evidence for 
substructure (thin and thick line, respectively). Bottom panels: 
same as top panels, but for M lT in lieu of M v . The masses M v 
and M, T are estimated from projected distributions of DM par- 
ticles. Initially, 500 DM particles are randomly selected along 
three orthogonal line-of-sights in each cluster, within a cylin- 
der of 1.5 /i _1 Mpc radius. Interlopers are then rejected using 
the procedure described in Sect. [3] Note that the vertical range 
is different from that of Fig.[3]and the vertical axis is now log- 
arithmic. 

when either a compact dynamical system is artificially split in 
two by the weighted gap algorithm, or, more frequently, a phys- 
ically distinct group is merged together with the main cluster. 
This last case is illustrated in Fig. [9] where we plot the ve- 
locity distributions of i) the DM particles within a sphere of 
1.5 h~ l Mpc radius (hatched histogram), ii) all DM particles 
within the cylinder of same radius (thin line), and iii) the DM 
particles identified as cluster members (thick line). In one of 
the three projections (bottom panel of Fig.|9j the cluster veloc- 
ity distribution is heavily contaminated by that of a foreground 
group, and the interloper rejection technique fails to isolate the 
cluster members. Note, however, that large errors in the mass 
estimates generally occur only for one of the three orthogonal 
projections considered per cluster. 



4.2. Subclustering 

The presence of subclustering has long been thought to poten- 
tially cause incorrect cluster mass estimates (e.g. Bird 1995 1. 



o 
o 
o 




0.4 0.6 0.8 1.0 
R (h- 1 Mpc) 

Fig. 7. Line-of-sight velocity vs. projected clustercentric dis- 
tance for a subset of the DM particles in the stacked clus- 
ter (see text). Filled squares indicate DM particles contained 
within a sphere of 1.5 h Mpc radius and selected as cluster 
members by the interloper rejection procedure, open squares 
indicate DM particles located within this same sphere yet in- 
correctly rejected as interlopers, and x symbols indicate DM 
particles located outside this same sphere and yet selected as 
cluster members (interlopers). 
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Fig. 8. The projected spatial distribution of the same DM parti- 
cles shown in Fig.0 Symbols are the same as in Fig. [71 



Several algorithms have been developed over the years to ad- 



dress this problem (see, e.g., Girardi & Biviano |2002 and 
references therein). One of the most commonly used (and 
easy to implement) algorithms is that of Dressier & Shectman 
( 1988; DS hereafter). Here we apply this algorithm (or, actu- 
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Fig. 9. The velocity distributions of a random subset of DM 
particles of a simulated cluster, along three orthogonal projec- 
tions. Thin lines indicate the velocity distributions of all DM 
particles within a cylinder of radius 1.5 h~ x Mpc, thick lines 
the velocity distributions of those DM particles identified as 
cluster members, and the hatched histograms show the velocity 
distributions of the DM particles located within a sphere of 1.5 
h~ l Mpc radius. 



Fig. 10. The ratio of estimated to true dynamical quantities, av- 
eraged over all cluster projections, as a function of the num- 
ber of DM particles selected as cluster members. Upper panel: 
M v /M v (lines above) and M a /M v (lines below). Lower panel: 
j r a, P lr a (lines above) and V3o- fl;P /cr fl (lines below). Dotted 
lines represent l-cr c.l. on the averages. The symbols represent 
the values obtained using galaxies (rather than DM particles) as 
tracers. Empty symbols are for all galaxies; filled symbols for 
galaxies with zj > 1.25 only. Stars refer to M v /M v and jRh/r/,, 
squares to Ma-/M v and <t 2 oo/o",,. 



ally, Biviano et al.'s 2002 version of it) to the simulated clus- 
ters seen along each of three orthogonal projections. As before 
(see Sect. 14. 1 i 500 DM particles are randomly selected within 
a cylinder of 1.5 A -1 Mpc radius; the DS test is applied on those 
particles that are selected as cluster members. 

In 59 cases out of 186 the DS test provides a probability 
> 0.99 that the cluster, as seen along the chosen projection, has 
significant subclustering. The cluster projections with > 0.99 
subclustering probability are also those with the largest mass 
overestimates. As an example, in the case-study illustrated in 
Fig. |9] subclustering is detected with high significance along 
the projection shown in the bottom panel. Among the 127 clus- 
ter projections without evidence for subclustering, only one 
has its virial mass overestimated by a factor > 2. As a con- 
sequence, the virial mass estimates are now closer to the true 
cluster masses, on average (see Fig.|6]and Tabled Id. no. 2). 

Hence the simple DS test appears to be a useful tool 
for identifying those cases of large virial mass overestimates. 
However, the situation is not so simple. When only 100 DM 
particles are selected instead of 500, the DS test gives a prob- 
ability > 0.99 only for 28 out of the 186 cluster projections. 
Clearly, when the sample size is reduced the DS statistics is less 
significant. Yet, observational results based on even smaller 
sample sizes, suggest a higher fraction of clusters with subclus- 
tering, ~ 1/3 (e.g. DS; Biviano et al. 1997 1. Hence the simu- 



lated clusters seem to display a lower amount of subclustering, 
on average, than real clusters. 

4.3. Undersampling 

In order to study how the accuracy of mass estimates depends 
on the number of cluster members available, we randomly 
select from 10 to 500 DM particles within cylinders of 1.5 
ft -1 Mpc radius, along three orthogonal projections. Cluster 
members are then selected using the interloper rejection pro- 
cedure described in Sect. [3] 

Reducing the sample size increases the scatter of the dy- 
namical estimates and affects the virial mass estimates (see 
Fig.[n|]and Tabled Id. nos. 1, 3, and 4). The average r ap jr a 
ratio tends to increase with decreasing number of selected clus- 
ter members, N m , for N m % 60, until it reaches a plateau for 
N,„ 55 20. On the other hand, the average V3cr fl p /cr fl ratio 
hardly changes with the size of the sample. As a consequence, 
for samples of ^ 60 cluster members, the simpler estimate 
appears to be a better predictor of the true mass than the virial 
mass estimate, M v , which is affected by the bias of r a ^ p . 

The increasing overestimate of the harmonic mean radius 
(and hence of the virial mass) as the size of the sample de- 
creases is mostly due to an increasing fraction of interlopers, 
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from 18% for N m £ 60, to 23-25% for N m £ 40. Part of the 
effect is however induced by a statistical bias that affects the 
harmonic mean radius estimate for small sample sizes. This 
can be seen by computing the harmonic mean radius in projec- 
tion, but using only real cluster members. We find that in this 
case, the ratio r a>p lr a remains close to unity for N m ^ 60, then 
starts increasing with decreasing N m , until it becomes ^1.1 for 
N m * 20. 

4.4. Radial-dependent incompleteness 

When a cluster field is observed in spectroscopy, complete- 
ness to a given limiting magnitude is rarely achieved. However, 
apart from the very bright cluster galaxies, there is no evidence 
for luminosity segregation in galaxy clusters (e.g. Biviano et 
al.lT992l 120011 Adami et al. I1998al Cappi et al. 120031 Goto 
2005 , Yang et al. 2005 ). Hence a random selection of the avail- 
able galaxies in a cluster field, even if not complete down to 
a given magnitude, should produce a sample with an unbi- 
ased projected phase-space distribution. In the previous anal- 
yses (Sects. |4~TI 14.21 l4~3t such a random selection was indeed 
applied. However, this is not always an easy task to accom- 
plish observationally. In fact, the number of slits (fibers) per 
mask in multi-slit (multi-fiber) spectrographs is fixed, hence 
the central, high-density regions of galaxy clusters are often 
sampled to a brighter magnitude than the external regions. As 
a consequence, the spatial distribution of the galaxies selected 
for spectroscopy turns out to be more extended than the parent 
spatial distribution of cluster galaxies. Said otherwise, a higher 
fraction of galaxies is selected from the total in the outer re- 
gions, i.e. the incompleteness is not random but depends on 
radius. 

Ignoring the problem of radial-dependent incompleteness 
can have quite a catastrophic effect on the cluster mass esti- 
mate, as we show in the following. We simulate an observa- 
tional set-up in which our clusters are observed with a multi- 
slit spectrograph with four quadrants, each 13 arcmin on a side, 
with a typical slit separation of 10 arcsec (we remind the reader 
that the clusters have been set at a distance of 250 h' 1 Mpc 
from the observer). The results are listed in Tabled Id. no. 5. 
In comparing these results to those obtained in the case of a 
random selection of tracers (Tabled Id. no.l), we note three 
differences. First, the average number of selected DM particles 
is only 198 out of the initial 500 randomly chosen in a cylin- 
der of 1.5 hT l Mpc radius. Second, the fraction of interlopers 
is increased (26% vs. 18%). Third, the harmonic mean radius 
is strongly overestimated in projection. All these differences 
are due to the geometrical constraints imposed by the multi-slit 
single-mask observation, that forces a sparser sampling of the 
dense cluster core relatively to the cluster outskirts. The obser- 
vational spatial distribution of the selected tracers is thus less 
centrally concentrated than the underlying parent distribution. 

In order to avoid this bias, the central cluster regions must 
be observed with more masks per unit area than the external 
cluster regions. Alternatively, photometric observations can be 
used to estimate, and correct for, the radial incompleteness of 
the spectroscopic sampling. When neither option is viable, it 



is better to rely on the M a estimate, which is much less af- 
fected by problems of incompleteness since cluster velocity 
dispersion profiles have only a mild radial dependence (see, 
e.g., Girardi et al. H996l and Fig.gJ. 

4.5. Other observational effects 

Other observational effects that could, in principle, affect a 
cluster mass determination, are i) a different size of the obser- 
vational aperture (the radius of the cylinder within which the 
tracers are selected), ii) the uncertainty on the determination 
of the cluster centre, and iii) errors on the redshifts of cluster 
galaxies. All these effects can change the projected phase-space 
distributions of the tracers, thereby affecting in principle also 
the identification of cluster members. 

In order to estimate the effect that different aperture sizes 
have on the cluster mass estimates, we compare the results ob- 
tained by selecting DM particles within cylinders of 1.5, 0.75, 
and 0.5 hr l Mpc radii. As the aperture a of the cylinder is de- 
creased, more particles are selected as cluster members out of 
the initial 500 distributed in the cylinder, but a higher fraction 
of the selected members are actually located outside a sphere of 
radius a. As a consequence, the results of the dynamical anal- 
ysis are similar for different apertures (compare the results for 
an aperture of 0.5 h~ l Mpc with those for an aperture of 1.5 
h Mpc, Id. nos. 6 and 1, respectively, in Table[Q. 

In order to simulate the observational uncertainty in the de- 
termination of a cluster centre, a random offset is added to the 
position of the real cluster centre (see Sect. [3}- The offset is 
randomly taken from a lognormal distribution with average of 
40 kpc and dispersion of 50 kpc, modeled after the observed 
distribution of uncertainties in the centre positions of ENACS 
clusters (Adami et al. 1998b 1. We find that centering errors of 
this size have essentially no effect on the accuracy of mass es- 
timation. 

Velocity errors contribute to increase the estimate of the 
cluster velocity dispersion. In order to simulate the observa- 
tional errors on galaxy redshifts, a random velocity offset is 
added to the real DM particle velocity. The offset is randomly 
taken from a gaussian distribution with width of 70 or 210 
km s -1 . An error of 70 km s _1 is typical for observations of 
nearby clusters observations (see, e.g., Katgert et al. [T996 1, 
while an error three times larger corresponds to the observa- 
tional uncertainties of galaxy redshifts in distant clusters (e.g. 
Demarco et al. 2005 1. As long as the velocity dispersion esti- 
mate is corrected following the prescriptions of Danese et al. 
(1980 1, the dynamical estimates are not significantly affected 
by errors on the velocities of the tracer particles. Note however 
that the 62 clusters used in our analysis have a median line-of- 
sight velocity dispersion of 622 km s _1 . Galaxy velocity errors 
are likely to become more important for mass estimates of low- 
velocity dispersion groups. 

4.6. Galaxies 

The limited number of simulated galaxies per cluster makes 
it impossible to study in detail all kinds of different observa- 
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Fig. 11. Same as Fig. [6]but for galaxies instead of DM parti- 
cles. No distinction is made here between projections with and 
without detection of substructures. 

tional effects, as it was done for the DM particles. But galax- 
ies have the advantage that early- and late-types can be distin- 
guished, based on their formation redshift z/ (see Sect.|2j. We 
start by considering all available galaxies within a cylinder of 
1.5 hT x Mpc radius. The results are shown in Tabled Id. no. 
7 and Fig.[^ Dynamical estimates based on galaxies are sig- 
nificantly affected by the presence of interlopers. About 1/4 of 
the selected cluster members are in fact outside the sphere of 
1.5 h~ x Mpc radius. As a consequence, both r a p and M v over- 
estimate the harmonic mean radius and the mass, respectively, 
by about 40%. In about 1/3 of the projections, the virial mass 
estimates are wrong by a factor of two or greater. On the other 
hand, the cr a p and M a estimates are (on average) almost correct 
(to within 5-6%). 

In order to compare these results with those obtained using 
DM particles as tracers of the potential, we randomly select 
within each cylinder of 1.5 h Mpc radius a number of DM 
particles identical to that of the galaxies in that same cylinder. 
In this way we make sure that the effects of undersampling (see 
Sect. 14. 3> are the same for the samples of galaxies and DM par- 
ticles. The results are given in Tabled Id. no. 8. The differences 
are very marginal. 

When only early-type galaxies (i.e. those with Zf > 1.25) 
are selected, the results of the dynamical analysis are quite dif- 
ferent from the case in which all galaxies are selected (compare 
Id. nos. 9 and 7 in Table ^respectively; see also Fig. XYl\ . In 
particular, the average fraction of interlopers is now reduced to 
18%, and the average value of M v is now almost identical to the 
real mass value. On the other hand, the M a estimate is 20% too 
low, as a consequence of a dynamical segregation effect (see 
the discussion below). The number of clusters with a mass esti- 
mate wrong by a factor of at least two is still high, about 1/3 of 



Fig. 12. Same as Fig.^^but only for galaxies with Zf > 1.25. 



the total, but this is expected because of the large scatter in the 
mass estimates due to the small number of available tracers. 

As before, we compare these results with those obtained by 
selecting the same number of DM particles as early-type galax- 
ies, within the cylinder of 1.5 \C l Mpc radius. The results are 
listed in Tabled Id- no. 10, and are quite different from those 
obtained using early-type galaxies, since, in particular, the frac- 
tion of interlopers is higher, and M v significantly overestimate 
the true mass. 

In order to understand the differences found between the 
mass estimates obtained using DM particles, all galaxies, and 
early-type galaxies, we consider their relative distributions in 
full phase-space. For this we build a stacked cluster, as in 
Sect. 13 The 3-d number density profiles of DM particles and 
galaxies in the stacked cluster are shown in Fig. ^] Galaxies 
are clearly less centrally concentrated than DM particles. This 
is confirmed by a Rank-Sum (RS hereafter) statistical test that 
reject the null hypothesis that the DM particles and the galaxies 
have the same radial distribution with a high confidence level 
(c.l. in the following), > 0.999 (see also Nagai & Kravtsov 
2005). Interestingly, this is in agreement with the observations 
that show a decreasing mass-to-light ratio as a function of clus- 
tercentric distance (e.g. Biviano & Girardi 2003 Rines et al. 
2004). On the other hand, the distribution of the normalised 
velocities, | v | /cr v , of galaxies is not significantly different 
from that of DM particles (RS test c.l. 0.829; see Fig .FBI. 

The spatial distribution of early-type galaxies is signifi- 
cantly more concentrated than that of all galaxies (RS-test c.l. 
> 0.999), and marginally more concentrated than that of DM 
particles (RS-test c.l. 0.939; see Fig.ll3>. The velocity distribu- 
tion of early-type galaxies is significantly different from those 
of all the galaxies and the DM particles (RS-test c.l. > 0.999; 
see Fig. 1141 . The high significance of this difference is a con- 
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Fig. 13. The 3-d number density profiles of DM particles (solid 
line), galaxies with Zf > 1.25 (dots), and all galaxies (squared 
x's) in the stacked cluster. 
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Fig. 14. The cumulative distributions of the 3-d velocities of 
DM particles (solid line), galaxies with zj > 1 .25 (dashed line), 
and all galaxies (dash-dotted line) in the stacked cluster. 



sequence of the large size of the data-sets, and does not reflect 
a large velocity bias: on average the modulus of the velocities 
of early-type galaxies is 0.94 that of DM particles and 0.93 that 
of all the galaxies. In real clusters, the different velocity distri- 
bution of early- and late-type galaxies is a well established ob- 
servational fact (Moss & Dickens fT9T71 Sodre et al. [T9l9l . The 
factor by which the velocity dispersion of early-type galaxies 
is lower than the velocity dispersion of all the cluster galaxies, 
is 0.93, on average, for the ENACS clusters (as can be derived 
from the estimates given in Biviano et al. 1997 1, a value iden- 
tical to the one found here for simulated galaxies. 

The different spatial distributions of all galaxies and DM 
particles does not seem to influence the results of the dynam- 
ical analysis in projected phase-space. One would expect the 
wider spatial distribution of galaxies to result in a larger r a , and 
hence r a p estimate. However, the r ap estimate is dominated by 
the presence of interlopers. Since the fractions of interlopers 



among galaxies and among DM particles are similar, also the 
r a p estimates are similar. 

When only early-type galaxies are selected, the fraction of 
interlopers decreases because the fraction of early-type galax- 
ies outside clusters is low, 29%, much lower than the fraction 
of early-type galaxies within clusters (see Sect. As a con- 
sequence, also the cluster r a p estimate is closer to the real r a 
value, than in the case of DM particles, although still ^ 20% 
too high. The r a overestimate does not result in an overestimate 
of the virial mass because cr ap underestimates the real cr a by 
^ 10%. This underestimate is caused by the narrower velocity 
distribution of early-type galaxies relatively to DM particles. 

In conclusion, galaxies have a biased distribution relatively 
to DM particles. Considering all the cluster galaxies, they have 
a wider spatial distribution than DM particles, but a similar 
velocity distribution. Considering instead only "early-type" 
galaxies, these have both a narrower spatial distribution, and 
a narrower velocity distribution than DM particles. Differences 
in spatial distributions have no effect on the dynamical mass 
estimates, since the estimate of the harmonic mean radius is 
dominated by the influence of interlopers. Selecting early-type 
galaxies helps improving the accuracy of the virial mass esti- 
mate because less interlopers enter the sample. On the other 
hand, the different velocity distribution of early-type galaxies 
with respect to that of DM particles, results in M a estimates 
that are too small. 



5. Discussion 

We have used a sample of 62 galaxy clusters extracted from a 
cosmological simulation, each with at least 10 galaxies within 
the cluster virial radius, to study the reliability of cluster mass 
estimates based on the distribution of their member galaxies. 
Galaxies in the simulation have been typed 'early' or 'late' 
based on the average formation redshift of their stellar popula- 
tion. The projected phase-space distributions of both the early - 
and the late-type galaxies of the simulation have been shown to 
be similar to those observed in real galaxy clusters (see Sect.|3. 

Two mass estimators have been considered in our analysis. 
One is the classical virial theorem estimate (corrected for the 
surface pressure term), and the other is an estimate based en- 
tirely on the cluster velocity dispersion (M, r , see Sect. [3}- By 
application to the set of simulated clusters in full phase-space, 
we have shown these estimators to be unbiased, and the clus- 
ters to be virialized (on average) as expected for z ~ clusters 
in a ACDM cosmology. 

In order to study how efficient are these mass estimators 
when applied to observed clusters, we have analysed our set 
of 62 simulated clusters in projection. Three orthogonal pro- 
jections have been considered for each cluster. DM particles or 
galaxies have been selected in cylinders of given aperture ra- 
dius and =i 200 h' 1 Mpc depth, in order to simulate the effect 
of interlopers. Cluster members have then been selected with 
methods commonly used in the recent literature (see Sect. [3}- 
DM particles have been considered as tracers of the poten- 
tial, instead of galaxies, when we wanted to explore the ef- 
fects of particular observational conditions with a sufficiently 
good statistics. However, the results of the dynamical analysis 
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are very similar when using the same number of DM particles 
or galaxies, except when galaxies of a given type are specifi- 
cally selected. Hence we think that our conclusions based on 
the analyses of clusters of DM particles should be applicable to 
the real world, when no distinction is made among galaxies of 
different types. 

Projection effects significantly affect the reliability of clus- 
ter mass estimates, through the inclusion of interlopers among 
the samples of presumed cluster members (see Sect. 14. 1> . They 
are responsible for large errors in the mass estimates of some 
clusters, depending on the projection direction. However, de- 
spite few catastrophic cases, the virial mass estimates are, on 
average, within 10% of the true values, if a sufficient num- 
ber of DM particles are used as tracers of the potential, k 
60. Nowadays, having a spectroscopic sample of <: 60 mem- 
ber galaxies per cluster does not represent a challenging ob- 
servational target, even for distant clusters (Demarco et al. 
2005 1, thanks to the multiplexing capabilities of instruments 
like VIMOS at the VL T (e.g . Czoske et al. l2002l and LRIS at 
Keck (e.g. Goto et al. 2005 ). This was not the case for large 
spectroscopic surveys of galaxy clusters of the past (Katgert et 
al. lT996l Carlberg et al. H997bi . 

Subclustering and the presence of groups along the line- 
of-sight to a cluster are responsible for the few catastrophic 
cases of strong mass under- and overestimate. Using the sim- 
ple DS-test for subclustering it is however possible to identify 
the worst cases and eventually remove them from the sample 
(see Sect. 14. 2> . Similar conclusions have already been reached 
by vH97 on their simulated clusters, and by Bird (1995 1 ob- 
servationally. Note however that it becomes increasingly dif- 
ficult to identify clusters whose masses are overestimated be- 
cause of subclustering, when the size of the data-sample is de- 
creased. Even with ~ 100 tracers of the gravitational potential 
per cluster, 3/4 of the cases of significant subclustering previ- 
ously found using ~ 500 tracers, can no longer be identified. 

The (positive) bias of virial mass estimates strongly in- 
creases as the number of tracers of the potential is decreased 
below ~ 60 (see Sect. 14. 3> . Since cluster velocity dispersions 
are always slightly underestimated (by ~ 5%) the bias has to 
do with an overestimate of the harmonic mean radius. This in 
turn has its origin in the presence of interlopers among the par- 
ticles identified as members. As already noticed by C97, these 
interlopers are characterized by a kinematic component of then- 
velocities that cancels the difference in the cosmological com- 
ponents of the particle and cluster velocities, i.e., these DM 
particles are infalling into the cluster. Since they are outside the 
cluster, their spatial distribution is not as centrally concentrated 
as that of real cluster member particles, hence they contribute 
to increase the harmonic mean radius estimate. Diaferio et al. 
dl999> came to essentially the same conclusions. 

In order to reduce the mass discrepancies when 55 60 mem- 
bers per clusters are available, one could take advantage of the 
fact that the estimates become less biased as the size of the 
sample is decreased. I.e., if the sampling of a cluster is poor, it 
is better to use rather than the virial mass. This is also the 
case when there is considerable uncertainty in the complete- 
ness of the selected sample of tracers (see Sect. 14. 4> . Sampling 
the same number of tracers within denser cluster regions (i.e. 



reducing the aperture size of the observational set-up) does not 
significantly improve the virial mass estimate, but could be ob- 
servationally more convenient. 

When using galaxies instead of DM particles as tracers, 
it is possible to draw subsamples selected on the basis of the 
galaxy properties. When only early-type galaxies are selected, 
the bias in the virial mass estimate is strongly suppressed, even 
if the average number of cluster members is very small. A sim- 
ilar conclusion was reached observationally by Biviano et al. 
J1997L who suggested to exclude emission-line galaxies from 
the sample of objects to be used in a cluster virial mass deter- 
mination (see also Sanchis et al. 2004 1. The improvement in 
the virial mass estimate obtained using only early-type galax- 
ies depends on the fact that the fraction of early-type galaxies 
is higher among cluster members than in the field, the so-called 
morphology-density relation (Dressier 1980) which is evident 
also in our simulated clusters. Hence, the interloper contami- 
nation is substantially reduced, compared to the case in which 
all galaxies (or DM particles) are selected. 

A comparison of our results with previous works is not 
straightforward, because of several differences in the analy- 
ses. Some of the previous works did not use the virial mass 
estimates at the estimated virial radius, but only the cr v esti- 
mate (Frenk et al. 11990; vH97), or isothermal mass estimates 
derived from cr v (C97; RB99), without using the spatial distri- 
bution of cluster members (i.e. the harmonic mean radius esti- 
mate). The M, T estimate used here differs from the isothermal 
mass estimates used in C97 and RB99 because their masses 
are derived within fixed linear radii, and hence have a crj de- 
pendence, while M, T is the mass within r v , and hence has a crj 
dependence. Both Sanchis et al. {2004 1 and Lokas et al. (2005 1 
derived the cluster masses from the application of the isotropic 
Jeans equation. Since this method requires a sufficient number 
of tracers, Sanchis et al. (2004i and Lokas et al. (|2005l only 
considered the case where 400, respectively 300, DM particles 
are available in each simulated cluster. 

Most previous analyses were based on the classical 3-cr 
clipping method of Yahil & Vidal J1977I . while the algorithm 
for the identification of interlopers used in this work is more so- 
phisticated. However, the methods used by Sanchis et al. ( 2004 1 
and Lokas et al. (|2005 1 have several similarities with the one 
used in this paper. Moreover, vH97 did consider the method of 
den Hartog & Katgert ( 1996 1, which is part of the technique 
used in this paper (see Sect, yj; vH97 indeed concluded that 
this method performs better than Yahil & Vidal's, particularly 
in the presence substructures. We do not expect to see a radical 
change in the average performance of the virial mass estima- 
tor as a consequence of using a different interloper rejection 
technique. First of all, it has been shown by several authors 
that different interloper rejection techniques generally lead to 
similar cluster mass estimates, when the number of sampled 
galaxies is sufficiently large (e.g. Girardi et al. 1993 Adami et 
al. U998d Biviano & Girardi 2003 Pimbblet et al. l2005l . Since 
a cluster velocity dispersion decreases with radius (see Fig.Q 
see also den Hartog & Katgert [19961 Biviano & Katgert[2004 1 
an interloper rejection method based on galaxy velocities only 
(like, e.g., Yahil & Vidal's) will probably reject more cluster 
members near the cluster centre, while accepting more inter- 
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lopers at the cluster edge, relative to a method that takes into 
account galaxy velocities and positions. However, the radial 
gradient of a cluster velocity dispersion is not strong, and the 
effect of neglecting the radial variation of a cluster velocity dis- 
persion on the estimate of global cluster quantities, like its mass 
and velocity dispersion, is likely to be negligible, unless there 
is substantial substructure. Substructure is definitely identified 
much better by using the combined spatial and velocity infor- 
mation than with velocity information alone (see, e.g., Girardi 
& Biviano 2002 and references therein). A detailed comparison 
of various methods of interloper rejection is however deferred 
to a forthcoming paper (Girardi et al., in preparation). 

Both Sanchis et al. (I2004t and Lokas et al. (2005 1 con- 
cluded that the virial masses of clusters can be reliably esti- 
mated using several hundreds of tracers of the gravitational po- 
tential, and we substantially agree on this. The method they 
used for estimating cluster masses are however different from 
the methods we considered in this paper. Both Sanchis et al. 
(20041 and Lokas et al.(2005 see their Figure 2) found that 
their method yields a mild systematic underestimate of the 
cluster masses. On the other hand, for similar numbers of trac- 
ers, our methods provide mass estimates within ^ ±10% of 
the true values. Their method is however aimed at deriving not 
only cluster masses, but also internal velocity anisotropics of 
the tracers population. 

The role of interlopers in the overestimate of the harmonic 
mean radius was emphasized by C97 and Diaferio et al. (Q999 1, 
and is in agreement with our finding. The systematic underes- 
timate of <t v by den Hartog & Katgert's (1996 1 method was 
already noted by vH97, but the effect was claimed to be larger 
than we see here, ~ 100 km s vs. our estimate of only 50 
km s _1 on average. C97 and vH97 suggested that the mass esti- 
mate bias decreases with the mass of the system, and increases 
with the aperture used to select the tracers. In our analysis we 
find that these trends are very marginal, if they are present at 
all. Perhaps this is due to our use of a more efficient interloper 
rejection algorithm. 

The issue of a velocity bias of galaxies relatively to DM 
has been raised several times, with some authors claiming it 
to be rather strong (Frenk et al. 19961, and others rather mild 
(Berlind et al. 2003 Faltenbacher et al. l2005> . We find that the 
result depends on which galaxies are selected, since the bias is 
negligible when all galaxies are considered, but significant, al- 
though small (0.94, on average), when only early-type galaxies 
are considered (see Sect. 14. 6> . The dependence of the velocity 
bias with galaxy mass will be explored in a future paper. 

6. Conclusions 

Using a set of 62 clusters extracted from a ACDM cosmolog- 
ical simulation, we investigated the behaviour of mass estima- 
tors based on the projected phase-space of galaxies. Several 
observational effects were considered, like the presence of in- 
terlopers and subclusters, the sample size, the size of the ob- 
servational aperture, incompleteness, and the selection of dif- 
ferent tracers of the gravitational potential (early-type galaxies, 
all galaxies, DM particles). 
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Our results show that the virial mass estimator is almost 
unbiased for samples of ^ 60 tracers. The average bias is 1.10, 
with a scatter of 0.30 (0.40) for samples of 400 (60) cluster 
members. For smaller data-sets, the bias of the virial mass esti- 
mator increases with decreasing sample sizes, reaching a max- 
imum of =! 1.5-1.6 when the number of cluster members is 
decreased to ^ 15-20. The virial mass estimates can be im- 
proved by removing clusters with significant evidence for sub- 
clustering, or by selecting early-type galaxies as tracers. The 
M lT estimator, based entirely on the cluster velocity dispersion, 
has a bias of ^ 0.90 for sample sizes ^ 60 members, and is es- 
sentially unbiased for smaller data-sets. The scatter of the Mo- 
estimator increases from ~ 0.30 for samples of ^ 60 cluster 
members, to ~ 0.60 for samples of only 10 cluster members. 
Since early-type galaxies have on average a lower velocity dis- 
persion than DM particles, the M a estimator based on early- 
type galaxies only, is biased low. 

Our analysis therefore suggests that the distribution of 
galaxies in projected phase-space can be used to provide re- 
liable estimates of a cluster mass. In order to optimize the mass 
estimate, early-type galaxies should be preferentially used as 
tracers, if a virial mass estimate is required. A simpler, more 
robust, and less biased estimate, especially for small data-sets, 
can however be obtained from an estimate of the velocity dis- 
persion of all galaxies identified as cluster members. However, 
it is not the scope of this paper to set up an optimized obser- 
vational strategy for the determination of cluster masses. The 
interloper rejection technique adopted in our analysis has been 
proven to work rather efficiently, but perhaps there is room for 
improvement. We are going to compare different interloper re- 
jection techniques in a forthcoming paper (Girardi et al., in 
preparation). Further progress is to be expected by increas- 
ing the resolution of the simulations to check the stability of 
the phase-space sampling of galaxies, and by a more realistic 
galaxy classification, e.g. based on colours. 
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